Star polymers: A study of the structural arrest in presence of attractive interactions 
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Simulations and Mode-Coupling Theory calculations, for a large range of the arm number / and 
packing fraction n have shown that the structural arrest and the dynamics of star polymers in a 
good solvent are extremely rich: the systems show a reentrant melting of the disordered glass nested 
between two stable fluid phases that strongly resemble the equilibrium phase diagram. Starting from 
a simple model potential we investigate the effect of the interplay between attractive interactions 
of different range and ultrasoft core repulsion, on the dynamics and on the occurrence of the ideal 
glass transition line. In the two cases considered so far, we observed some significant differences 
with respect to the purely repulsive pair interaction. We also discuss the interplay between equilib- 
rium and non equilibrium phase behavior. The accuracy of the theoretical tools we utilized in our 
investigation has been checked by comparing the results with molecular dynamics simulations. 
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I. INTRODUCTION 

When a substance dissolves into another to form a true solution and the units dispersed through the solvent are 
much larger in size than the solvent molecules, we name it a colloidal dispersion. There are different classes of 
colloidal suspensions: molecules individually larger than lnmQ e.g. proteins, polysaccharides, polymers, as well as 
dispersions that arise when a number of small molecules associate together to form an aggregate, it is known, for 
example, that amphophilic molecules in a suitable solvent, above a critical concentration, can aggregate in micelles 
Varying the interactions between the mesoscopic constituent particles in colloidal dispersions results in a broad 
range of equilibrium and non-equilibrium fluid behavior. 

One strategy to study such systems involves summing over the solvent degrees of freedom, leaving an effective in- 
teraction between the center of mass of the macroparticles The resulting effective potential V(r) depends, in a 
complicated fashion, on the interparticle separation and it is in general state dependent 0| , @ • A special example in 
this context is given by the pair potential describing a star polymer solution. 

Star polymers (SP) can be considered as a generalization of polymer-coated colloidal particles in the limit where 
the number of monomers (N) per chain is large and the size of the central core is small with respect to the whole star 
extension. The repulsive interaction between star polymers at short distances increases very slowly as the interparticle 
separation decreases, namely in a logarithmic way. The number of chains chemically linked to a core influences the 
'softness' of V(r): stars with small arm number / may interpenetrate widely; in the limit case / = 1,2 the star 
polymer reduces to a simple polymer chain. Stars with very large / emphasize their colloidal nature, although only 
in the limit / — > oo the potential resembles the hard-sphere one. 

Recently solutions of star polymers have received attention in relation to several medical and industrial applications 
0. Moreover, in the last ten years, advances in macromolecular chemistry, leading to the synthesis of regular star 
polymers, have made it possible to explore the physics of very different model systems monodisperse in iVand/ Q, 0- 
Finally these systems are very interesting from a theoretical point of view in relation to their polymer-colloid hybrid 
character. Indeed the efficient design of new mesoscopic materials with properties intermediate between different 
classes of colloidal systems is a very important challenge in soft condensed matter physics. 

In the last decades the equilibrium and non-equilibrium phase behavior of SP solutions have been widely investi- 
gated, both theoretically H H EJ 

and experimentally ||| fl3L ITiL [Tlj . These system are an interesting example 
of a complex fluid for which the phase diagram has special features arising from the ultrasoft nature of the repulsive 
interaction: e.g. there exists a cutoff value of the functionality / below which the system is fluid for all densities 
(f c = 34) and for / > f c the phase diagram exhibits several unusual solid lattices as well as reentrant melting [TT| : 
SANS and SAXS experiments on solutions of many-armed stars above f c , e.g. ^||, have revealed different macro- 
crystal structures as we increase the density. The functionality-dependent bec- and fcc-solids 0, 0] , as well as the 
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reentrant melting transition |19| . have been experimentally observed in solutions of star-like block copolymer micelles. 
A second freezing transition observed in the same experiments can be interpreted as the freezing in a bco crystal |20j | . 

The dynamical properties of star polymers have been extensively investigated. Several studies, focused on star- 
polymers or star-like systems in athermal solvent with different arm numbers, have shown that it is quite difficult to 
nucleate a crystal: in many cases, mainly at high functionalities, the solutions display a gelation transition |2lLl22ll23 . 
l^iEHEBl- Molecular Dynamic (MD) and Brownian Dynamic (BD) simulation data for a large range of / and packing 
fraction r\ = ^p<? 3 show that the dynamics of star polymer solutions is extremely rich. In particular the ideal glass 
transition line, obtained by means of the Mode-Coupling Theory (MCT) |27|. displays a non-monotonic behavior as a 
function of 77 and /. This behavior has been connected to the rj and /-dependence of the effective hard core diameter 
of an equivalent HS system. The detailed comparison between theoretical predictions and simulations confirms the 
validity of the MCT approach to study the disordered arrested states in soft matter like colloids [2^, and in 
ultrasoft systems like star polymer solutions |30l |. In particular it has been confirmed that the modified hypernetted 
chain integral equation (MHNC) is a very good approximation to study static correlations in systems described by 
ultrasoft interactions not only in equilibrium |3l| but also in metastable states |30j . 

Recently a model potential has been suggested to describe SP solutions where, in addition to the excluded volume 
effects, attraction emerges due to dispersion or depletion forces. For this model the fluid-fluid phase diagram has been 
determined I3II using mean field theory and two fluid-state-theories, MHNC and the hierarchical reference theory 
(HRT) |32, l33t l34j for different /. This analysis shows that when the strength of the interaction is strong enough a 
fluid-fluid phase transition appears but the density-temperature coexistence curve bifurcates at a triple point into two 
lines of coexistence terminating at two critical points. This peculiar phase behavior is related to the unusual form of 
the repulsive contribution. 

Moreover it has been shown that self-organized structures, resulting from telechelic linear homopolymcrs and copoly- 
mers, similar to star polymers, can bridge by producing an effective attractive interaction leading to reversible ag- 
gregation of macromolecules (35|. These micelles are constituted by telechelic associative polymers which have the 
associating groups at the chain ends. Above a critical concentration the end groups associate in multiplets, forming 
flower-like polymeric micelles. At higher concentrations the process of bridging can lead to the formation of a transient 
gel or also induce macroscopic phase separation |36l 137^ . 

The aim of this paper is to investigate the effect of attractive interactions on the dynamics and on the occurrence of 
the ideal glass transition line, trying to emphasize the special features of the phase behavior arising from the ultrasoft 
nature of this repulsive effective interaction. We also qualitatively discuss the interplay between equilibrium and non 
equilibrium phase behavior. In this work different theoretical and numerical methods have been utilized (mean field 
theory, fluid state theories, mode-coupling theory, molecular dynamic simulations). 

We recall that in the last years a great number of studies focused on the dynamical behavior of short ranged attractive 
systems, which are characterized by a strong repulsive core besides the attraction. In particular when the range of the 
attraction becomes much shorter than the typical diameter of the colloids, phenomena like a reentrant glass transition 
or the existence of two different glassy phases emerge |3^, E^llicl. This peculiar behavior has been confirmed by a 
large number of simulations 0, l42l 1431 Plij and experiments Notice that, historically, these new findings 

have been predicted for the first time within MCT calculation and only on a second stage confirmed by experiments. 
Consequently it is clear that this kind of approach can be extremely useful also for different interaction models. It is 
interesting now to focus the attention on ultrasoft repulsion (typical for example of a SP solution) and an attraction 
with different range (which could be typical, for example, of depletion interactions) and to investigate the possibility 
of new features. 

This paper is organized as follows: in section[H]we present the general framework for our research. In particular, in 
subsection 111 Al we introduce the interaction model we chose to study star polymers in presence of attraction. Then in 
subsections iHBl HTCl and iHDl we describe the theoretical and numerical tools we used to study the structural prop- 
erties, the structural arrest and the diffusivity in dense star polymer solutions: MHNC, MCT and MD, respectively. 
In particular we discuss the application of the modified hypernetted chain integral equation to study the structural 
properties of star polymer solutions and some test of its accuracy when attractive interactions are taken into account 
in addition to the entropic contribution. 

In section llTTl we consider the effect of attraction on the slow dynamics and structural arrest of star polymer solutions. 
We carried out a mode-coupling theory analysis which allows us to locate the nonergodicity transition curve of the 
system, using as input the information on the structure obtained by MHNC. Our aim is to complete the picture of the 
phase diagram of a star polymer solution in the presence of attractive forces, investigating the d yna mics, for the values 
of the parameters which govern the intensity of the attraction extensively discussed in Rcf. |31|. Then we modify 
these parameters in such a way to consider attractive forces of shorter range, and we focus on the effect of these 
modifications on the properties of the glass state. In order to test this difference, molecular dynamics simulations 
have been performed and the diffusivity of SP fluid as been evaluated up to crystallization. In this section we again 
discuss the structure of the system very close to the glass phase. 
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Finally in section HVI wc discuss and summarize our results and we draw our conclusions. 

II. INTRODUCTION TO THE MODEL AND METHODS OF STUDY 

A. Effective pair interaction 

The effective pair interaction between star polymers with / arms in a good solvent is purely repulsive and for / > 10 
it reads as follows: 
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a is the corona diameter of the star and depends on the number of monomers of a single arm [48| . Witten and 
Pincus determined by scaling theory the explicit form of such interaction at short distances [ij; it has been also 
shown that the good agreement between theory and experiments significantly improves if the model interaction has, 
in addition to a logarithmic repulsive core, a long range interaction of Yukawa form E3j E3> This model 
interacti on g ives a good description of small angle neutron scattering (SANS) results on concentrated star polymer 
samples [l5| for an wide range of / values. For / < 10 a logarithmic-Gaussian potential more accurately describes 
the effective interaction |5J but we do not study this regime. 

In a previous paper |3lj the phase diagram of star polymer solutions has been investigated when the effective 
interaction between two star polymers contains an additional attractive contribution w(r) to be added to the repulsive 
part V rep (r) of equations 

V tot (r)=V rep (r)+w(r). (2) 

This attraction could stem, for example, from a van der Waals interaction arising from a non perfect matching between 
the refraction index of the solvent and of the polymer in a way that does not alter the basic configuration of the single 
star polymer; in this case w(r) is independent from the temperature. Alternatively attraction could be induced by 
depletion interaction when a third component, which is large compared with the solvent molecules but small compared 
with the star polymers, is present in the solution. In this case w(r) has an entropic origin and consequently it depends 
on T. In this situation the size and concentration of the depletant are, respectively, closely tied to the range and 
strength of the attraction. In what follows we shall use temperature as control parameter. 

To study the structural arrest and the dynamics in presence of attractive forces we can use the simplified model 
potential utilized in ref. |3l| where w(r) has the functional form of a Fermi distribution, i.e. 

C 

W ^ = lr-Al I -I " ( 3 ) 

exp [^\ + 1 

The parameters A and B control the position and the width of the well potential, C is the amplitude of the attractive 
contribution. For convenience we use reduced units for temperature in terms of C: T* = ksT/C. By a suitable 
choice of these parameters one can guarantee, for the temperatures of interest, that Vt i(r) does not have a significant 
subsidiary maximum at large r and avoid possible complications due to the competition over different length scales. 
We have studied three sets of parameters A and B in order to see the dependence of the glass transition on the 
range and the width of the attractive well. The first case corresponds to the potential discussed in Ref. [^H , i-e. 
A = 2.1a, B = 0.35cr and we refer to this case as SPi. For these values of the parameters, the attractive contribution 
is rather long ranged and we found that w(r) does not modify the ideal glass transition with respect to the purely 
repulsive case. For the other two choices of the parameter, SP2 corresponding to A = 1.95a and B = 0.21er and SP3 
corresponding to A = 1.875er and B = 0.155<r the attractive well is displaced to smaller distance and we observed 
significant differences for the glass transition with respect to the only repulsive pair interaction. In Fig. ^ we show 
the shape of the potential in the three cases for / = 50. It should be noticed that due to the ultrasoft character of 
V rep the value of A and B modifies not only the range of the interaction and the position of the minimum but also 
the depth of the well potential. In addition, the position and width of the attractive well depends on temperature 
because V rep scales with T. In Fig.^wc plot also V rep : at short distances we can observe that Vtot is softer than V rep . 
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B. Modified hypernetted chain integral equation 

In the present work we carry out a Mode-Coupling Theory (MCT) analysis [27| of the long-time limit of the 
correlation functions to locate the ideal glass transition line of the system. MCT provides a set of closed equations to 
calculate the non-ergodicity parameter f q which acts as an order parameter for the glass (sec III D|) . All the information 
needed to solve these equations is contained in the static structure factor, defined as S(q) = (p(—q)p(q)) /N and in 
the number density, p = N/V, p(q) being the density fluctuation variable of wave vector q. To calculate S(q), we 
utilized the MHNC integral equation. This equation is in general accurate also when an attractive contribution to the 
interaction is present [52] ■ In the case of star polymers in a good solvent it has been verified the remarkable accuracy 
of this theory, to describe fluid states as well as metastable states |3(| [3ll for a large range of / and density values. 

The starting point for the MHNC equation is an an exact relation |53l | . obtained from a cluster expansion, which 
connects the radial distribution function (rdf) g(r) to the interparticle potential V(r): 

g{r) = ejrp[-/3V(r) + h(r) - c(r) + E(r)] (4) 

where h(r) = g(r) — 1 and c(r) are the pair and the direct correlation function, respectively. c(r) is related to h(r) 
by the Ornstein-Zcrnike equation. 

The term E(r), called bridge function, represents a sum of an infinite number of terms, the so called elementary 
graphs in the diagrammatic analysis of the two-point function and, in general, it is not known. In the MHNC scheme 
the E(r) is replaced by the bridge function of a fluid of hard spheres, Eh sir) , of suitable diameter d. To optimize 
this choice, which depends on the parameter d, the free energy is minimized [54| . This is equivalent to satisfy the 
relation: 

[dr[ 9 (r)-g HS (r,r ms) ] ^(^hs) =Q } 
J arjHs 

where t]hs = %P<P 

In order to implement MHNC one needs the rdf <?irs(r) of hard spheres from which one can obtain Ehs- Vcrlet and 
Weis (VW) provided an accurate parametrization of gHs( r ) based on the PY equation with a correction which 
incorporates thermodynamical consistency through the Carnahan-Starling state equation |53l ] . This together with 
equation il llol) gives a closed set of equations which are solved by a standard iterative method. 

In a previous paper some of the authors showed that the dependence of tjhs on the density as determined by Eq. © 
is unusual and reflects the features of the interparticles interaction j^H . Moreover in ref . we found evidence that 
the characteristic sequence of maxima and minima of tjhs versus r\ is directly related to the non-monotonic behavior 
of the diffusion coefficient as a function of the packing fraction. Indeed the slow dynamics in star polymer systems 
at low and intermediate densities can be qualitatively described as the slow dynamics of the hard sphere system via 
a density and functionality dependent effective diameter determined with MHNC closure. 

C. Molecular Dynamic simulation: comparison between MHNC and MD results on the correlations. 

To check the results of the MHNC and MCT calculations, we performed extensive MD simulations for the model 
described by the potential defined by Eq. (lll.'il) . We simulate N = 1000 particles for different values of the density, 
temperature and functionality. For each state point the configurations have been equilibrated at constant temperature 
for a time long enough to ensure both the equilibration and decorrelation from the initial configuration. The acquisition 
run started after this preliminary preparation, then the system was simulated at constant energy. We carefully checked 
whether crystallization occurred or not during the run by inspection of the static structure factor, with the same 
modulus, not averaged over different directions of g-vectors. In fact to improve the averages in the calculations of 
S(q) at a given q, we generally considered up to 300 independent g-vectors chosen with a random direction but with 
the same modulus q. In correspondence of these values, we evaluated the density variables and consequently the 
static structure factor for a given direction. Finally we averaged on the different directions. For a liquid the structure 
is extremely disordered and all the terms for different directions give roughly the same contribution to the average. 
When the structure factor corresponding to a certain direction starts to grow over the others, the system begins to 
show a preferable direction, i.e. crystallization takes place. When a configuration crystallized we discarded it. As 
discussed for the purely repulsive case |30j . the system has a strong tendency to crystallize, being monodisperse in 
diameter [44| . Consequently, we will test our predictions only when the system remains in the liquid phase, a future 
direction of our research may be directed toward the suppression of crystallization introducing, for example, a second 
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component slightly different in diameter [4^ . We chose as unit of length the corona diameter a and as unit of mass 
the mass of the particles. Moreover we measure the temperature in reduced units, i.e. T* = kgT/C. 

We begin our discussion on the numerical simulation by discussing the purely repulsive case, i.e. C = 0. As 
extensively described in ref . [3lj , MHNC gives a very accurate description of correlations for star polymers in a good 
solvent for a wide range of densities. In Fig.[3]we present results for the repulsive case, for / = 32 and r\ = 0.6 obtained 
from MHNC and Rogers- Young (RY) theoretical calculations and from MC and MD simulations [^3 The RY equation 
is another integral equation for g(r) which interpolates between PY and HNC equations |57j . The agreement between 
MD and MHNC, and between MD simulations and MC simulations, is extremely good. The small apparent MD 
underestimation of the main peak is due to the grid we used in the simulations. A more refine grid, however, would 
increase the noise in the data. On the other hand RY shows some significant discrepancies with respect to simulations 
results in the range of strong coupling which emerges for packing fractions of order of rj = 0.60 and r\ = 3 |3l| . By 
comparison with simulation results (both MC and MD) we definitively conclude that MHNC describe accurately the 
structure of star polymers in a good solvent. 

Hence, we focused on the attractive case and we compare the MD results and the prediction of the MHNC theory 
for the fluid case: in particular we compared the radial distribution function and the structure factor for the three 
different models of the attraction described in section Til Bl in a regime of strong coupling when the main maximum of 
S(q) reaches large values. The temperatures and packing fractions we have chosen in our study, correspond to states 
of interest in the study of the ideal glass transition line, i.e. close to the fluid-glass and glass-fluid ideal transition line. 
We find in general a very good agreement between MD results and MHNC results for all the values of temperature, 
packing fraction and arm number investigated. In particular a good accuracy is achieved in the determination of the 
position of the main peak in the radial distribution functions and structure factors. Only for states at low temperature 
and at packing fraction corresponding to very large coupling (e.g. see Fig. |3J) one notices some small discrepancies 
between MHNC and MD results. 

This analysis confirms the accuracy of the MHNC integral equation in describing the structure of an ultrasoftcore 
potential also in presence of an attractive contribution at longer range. 

One can notice that the MHNC g(r) has some structure in the region of the second maximum and this can even 
show up as a subsidiary maximum (Fig. [3J. It is known |52l| that this spurious structure is a consequence of the VW 
parameterization of gns m terms of the PY solution but it is believed that this anomaly has no serious consequence 
in S(q). 



D. Mode Coupling Theory for the Ideal glass transition 



In this section, we shall briefly review the nature of MCT, and discuss the type of information it yields. The MCT 
of super-cooled liquids describes the dynamical transition by a nonlinear integro-differential system of equations for 
the normalized time correlation functions of density fluctuations <f>(g, t) = (p(—q,t)p(q,0)}/(p(—q 7 0)p(q,0)), where 
p(q, t) = J2i ex P (*Q ' r i(t))- As discussed above, the only input to the MCT equations are the equilibrium static 
structure factor, S(q) and the number density, p. The glass transition can be identified by studying the long time 
limit of the MCT equations, which determine the non-ergodicity parameter of the system f q = lim^oo $>(q,t). An 
crgodic state is characterized by f q = 0. This value is always a solution of the MCT long-time limit equations (271 ]. 

The quantity f q obeys the equation f q /(l — f q ) = J- q (f). Here, the mode-coupling functional T q is given by 



(27r) 3 q'kJkJ \q-k\ ■ 



(6) 



Equation [|| together with the equation for f q can be derived by taking the long time limit, t —* oo, of the MCT 
equations |27| . The mode-coupling vertices are determined by the structure factor S q , the direct correlation function 
c q , and the density p: 



V q,k = S l S kS w _i\P q-kc k + q- {q~k)c w _^ /q 



(7) 



The direct correlation function is directly related to the static structure factor by the relation c q = ^ [1 — -g-J. The 

glass transition appears as an ergodic to non-ergodic transition for the system, where f q y^0 solutions arise. These 
transition points correspond to bifurcation singularities of the MCT Eq. (|6I7|I . 

In the present work, we numerically solved Eq. I|6I7|I with an iterative procedure over a grid of 650 cqui-spaced q- 
vectors up to q = 62.32. For the static structure factor we used the static MHNC structure factor S(q) calculated as 
described above. 
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III. RESULTS 

A. SPi: MCT and MD results 

First of all we considered the model potential we labeled SP\. At low temperature the system has fluid- fluid 
phase transition with two critical points: the first critical point is around rj c = 0.026 and T* ~ 0.7, the second one 
around r] c = f.04 and T* ~ 0.f56. Far from the critical point, we observed that the structure factor S(q), relative 
to the potential V tot , is not very different from the structure factor relative to the simply repulsive interaction V rep . 
In particular the position of the main peak does not change appreciably and its height change for less than 1% 
(0.2 < T* < 0.6). Consequently one might expect that the location of the ideal glass transition will not be different 
from the repulsive case [3(j|. Indeed this is the case. We studied several arm numbers (from low value, / = 24,32 
to higher one / = 70) for several values of the packing fraction with MCT. Similarly to the case of purely repulsive 
interaction we deduce that below / = 46 there is no glass phase. For 46 < / < 60 the system is fluid for low densities 
and glass for intermediate (0.25 < rj < 0.7), but it is still fluid at high densities (77 > 0.7). When 60 < / < 72 the 
system is in a glass state for intermediate and very high densities (0.25 < 77 < 0.7, r\ > 2.25). Finally above / ~ 72 
the system is fluid for low densities (77 < 0.25) and in a glass state for intermediate and high densities (77 > 0.25). 

Hence, the effect of the attraction does not change the location of the ideal glass transition. Consequently, in the 
temperature-packing fraction plane, the ideal glass transition line would be trivially represented by a vertical line. 

We also performed molecular dynamic simulations to calculate the dependence of the diffusion coefficient, from the 
long time limit of the mean squared displacement, varying the temperature, for different values of / and the packing 
fraction. This analysis supported the MCT predictions, i.e. the diffusion coefficient does not change appreciably with 
the temperature. 

Starting from the behavior of the structure factor changing the temperature and from the analysis of the structural 
arrest we performed, we can complete the picture of the phase diagram. For / < 34 there is no freezing transition: 
the two fluid-fluid phase transitions represent stable states; for / > 50 all the density region above the density 
corresponding to the triple point T p is occupied by crystalline phases; for 34 < / < 50 bot the fluid-fluid phase 
transitions and their critical points persist as stable states. In Fig. 0] we describe the whole phase diagram for / = 50: 
the fluid-fluid phase diagram is calculated by means of mean field theory and the ideal glass transition lines determined 
with MCT. At the equilibrium the region around the second critical point becomes metastable with respect to the 
freezing, the regions between two successive squares (from low to high densities), in the upper part of the figure, 
indicate the densities where the system is solid (data are MC results reproduced by Watzlawek et al. ^3). If the 
crystallization can be avoided we can observe that the second critical point lies outside the glassy region (see also 
the inset). The limit of very high temperatures, i.e. the purely repulsive interaction are also shown in Fig. 0] the 
densities which delimit the glass region are effectively the same with and without attractive interaction. 

We conclude that, varying /, the first critical point will always lie outside the glass region. The second critical 
point, however, will enter the glassy phase for / > 70. 

B. SP 2 : MCT and MD results 

As anticipated above we tried to tune the parameters A and B in order to enhance the effect of the attraction on 
the shape of the glass-transition line. Our aim was to shrink the well and move the minimum of the pair interaction 
close to sigma. 

First of all we consider the parameters A = 1.95a an B = 0.21cr (SP2). As we see in Fig. Q] the well potential 
in this case is shrunken approximately by 30% with respect to SP± and the position of the minimum is changed 
approximatively by 6% (0.25 < T* < 0.6). If compared to SP\ the depth of the well potential changed from 5% with 
respect to at lower temperature till approximatively 28% around T* = 0.6. 

In Fig. [S] we present the occurrence of the glassy phase for different numbers of arm and as a function of temperature. 
For / > 46 we can summarize our results as follows: for high temperatures the two ideal-glass lines (fluid-glass on 
the left and glass-fluid on the right) tend to the repulsive case |30j . On lowering the temperature, however, a fluid- 
stabilizing effect sets in, so that the fluid-glass line tend to move to larger values of the density. So in the low density 
region, roughly below 77 = 0.50, there is the possibility of the glass melting when lowering the temperature. In this case 
increasing the temperature, the width of the nonergodicity parameter f q , which is a measure of the inverse of the cage 
localization length, gets larger (see Fig. HJJ). Moreover we observe an increase of f(q = 0) lowering the temperature, 
corresponding to densities closer the coexistence curve. The effect of the attraction for large value of / is very small 
and in particular for / = 70 we can see that the fluid-glass line is again very similar to the repulsive case. On the 
right side of the glass region we can observe that the glass-fluid line move to higher densities when the temperature 
decreases and consequently when the intensity of the attraction increases. In contrast to the low density case, the 
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attraction now favors the formation of the glass to higher density with respect to the repulsive case. On lowering the 
temperature the nonergodicity parameter f q presents a larger width in q. This clearly indicates that the particles are 
localized on a shorter length scale. This phenomenon could be related either to the effect of a stronger attraction or 
to a larger value of the effective repulsive length. It is interesting to stress that in all the cases taken into account 
the non-ergodicity parameter, presents the typical shape of a repulsive glass, i.e. a glass that possesses a structure 
dominated by the caging effect. This typical repulsive behavior is characterized by oscillations in correspondence 
to the peaks of the static structure factor with a maximum at the first peak. It has been shown, both by MCT 
calculations jijj and computer simulations 01 1 that for sufficiently short-ranged attraction a new glass, that has 
been named attractive glass, emerges. Indeed the shape of the f q for an attractive glass is completely different: the 
oscillations are very weak and the maximum is not so pronounced. In our investigation we encountered fingerprints 
only of the repulsive glass. The effect of the attraction and the interplay between attraction and repulsion, seem to 
act more on the size of the cages and on their formation rather than on the nature of the arrest itself. However we 
have not investigated the case of very narrow attractive wells on the length scale of a. 

The behavior of star polymers with / = 45 is different: we can observe that the attractive contribution of the 
interaction favors the occurrence of glass states at low temperature. Indeed for this value of the functionality, a stable 
glass phase emerges at a temperature T* < 0.6 from the solution of MCT equations and the density range of the 
glass phase gets larger on further decrease of the temperature. In this case the q-width of f(q) increases lowering 
the temperature, both for rj < 0.5 and ?y > 0.5. It is interesting to note that no glass transition is found in the 
corresponding repulsive case, i.e. in the limit of high temperature. 

In the so called "repulsive glasses" the occurrence of the ideal glass transition in the framework of the MCT depends 
strongly on the behavior of the main peak of the structure factor, i.e. on the first neighbor interactions. 
In Fig. we can observe a magnification of the main peak of the structure factor for / = 50 and two different packing 
fractions on the left and on the right side of the ideal glass region (MD simulations). For the different temperature 
investigated we observe the same trend: for r\ = 0.314 (Fig. 0-a)) decreasing the temperature the main peak of the 
structure factor decreases, showing a loss of the correlation between particles. The opposite trend is observed on the 
right side, r\ = 0.628 (Fig. 0-6)). We also studied the behavior of the first peak in the structure factor for / = 45, in 
this case the trend in the glass region is represented by an increase of the main peak lowering the temperature, indeed 
this effect is responsible for the anomalous formation of the glass at this value of the functionality. 

This preliminary and qualitative study, has been supplemented by a more thorough analysis: we performed molec- 
ular dynamic simulation for several values of /, T* and rj. 

In Fig. [S] we show the diffusion coefficient D/Dq as a function of the temperature calculated for packing fraction not 
far from the glass region and for different values of / (42, 45, 46, 50). In the present context data for D are normalized 
by Do = o\jT* /m, in order to take into account the temperature dependence of the microscopic time |44j . The 
MD-diffusion coefficients as a function of T* confirm the MCT trend. On the low density side of the glass region 
the diffusion coefficient decreases increasing the temperature while on the high density side the diffusion coefficient 
increases increasing the temperature. The values of D/Dq, for fixed temperature, decrease increasing /; we remark 
that for / = 42 the system does not show glass transition. 

In general it has been observed that the effect of the presence of a glass transition can be noticed as a decrease in the 
diffusivity also far from that part of the phase diagram where the structural relaxation time starts to grow. Indeed 
this has been encountered, for example, in a monodisperse square well system |44j and in a purely repulsive soft 
potential [3(j. In MD simulations both these systems presented a strong tendency to crystallize. However, studying 
the behavior of the diffusivity, it is possible to pin down the shape of the glass transition line. Also in the present 
work we managed to check that the diffusion follows the MCT predictions. The fact that the / = 45 trend could not 
be checked, i.e. an increase of the diffusivity on rising the temperature for both high and low density region, could 
be related to the fact that the simulations should be run closer to the glass transition and, unfortunately, this is not 
feasible due to the occurrence of the crystallization. One way around this problem would be to extend the results 
to a system that presents a smaller crystal nucleation rate. For the square well system, for example, this has been 
accomplished considering a binary mixture |42| . 

Summarizing we remark, for the different arm numbers investigated, that for ?y < 0.5 the attraction seems to 
determine a destabilization of the cages while for r] > 0.5 the attraction facilitate the formation of the glass for higher 
density with respect to V rep . The most peculiar behavior has been found for a value of / (/ = 45) for which the glass 
lines are very near to the value of 7/ — 0.5. For this packing fraction star polymers start to interpenetrate widely. 
We emphasize that, considering stars in a good solvent at the equilibrium, this value is just in the middle of the bec 
(34 < / < 55) or fee (/ > 55) crystal phase. Related to the behavior of the structure factor corresponding to SP2, 
we expect that also the densities where the system is solid, at the equilibrium, could change weakly with respect to 
the repulsive case. We will return more in detail on this point in the section relative to conclusions and discussion. 

Finally we show in Fig. |5] the position of the ideal-glass lines with respect to the fluid-fluid coexistence curves for 
/ = 50. It turns out that the temperature should be decreased a lot to notice possible effects of the density fluctuations 
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around the second critical point on the ideal glass transition line. 

C. SP 3 : MCT and MD results 

The third case we considered is characterized by A = 1.875(7 and B = 0.155<7, and we indicated it SP3. In this 
case the width of the well potential decreased, with respect to SP\, from approximatively 40% at T* = 0.25 till 60% 
at T* = 0.6. The position of the minimum is changed by less than 9%. If compared to SP\ the depth of the well 
potential changed from 7% with respect to at lower temperature till approximatively 48% around T* = 0.6. 

In Fig. EH the occurrence of the glass phase is shown for / = 46, 50, 70 for a wide range of temperatures. We 
analyzed in more detail the case / < 70 as here we expect very different behaviors from the three cases in exam. 

For / > 50 the qualitative behavior is very similar to that of SP2 ■ For densities below the glass region the influence 
of the attractive term, when lowering the temperature, is a liquid-stabilizing effect. In the opposite regime (densities 
above the glass region), increasing the intensity of the attractive contribution results in moving the glass- fluid line to 
higher densities with respect to V rep . We point out that here the curvature of the lines is more pronounced than in 
the SP2 case. Anyway comparing all the cases investigated is not trivial: this is in part due to the fact that a correct 
rescaling of the temperature should be done considering the change in the depth and the position of the well potential 
changing A, B and T, rather than the amplitude C of the attractive term. We will return on this point in section HVI 
For / < 46 the behavior of the star polymer solution is completely different from the previous cases. Indeed for / = 46 
it exists a value of the temperature, T* ~ 0.6, below which the system does not present a glass transition. For / = 45 
we do not observe any glass phase for all the temperature values investigated (0.4 > T* > 0.8). 

As for the previous case, we performed MD simulations for different arm numbers and temperatures, confirming 
the MCT result: the trend of D/D is qualitatively the same as for SP2 (see Fig. lTT|) . For / = 46 we do not capture 
the increasing of the diffusivity when decreasing the temperature, as we expect on the basis of the MCT results on the 
right side of the glass region. As in the previous case this might be due to the fact that we did not consider packing 
fraction close enough to the ideal glass lines. 

Finally we studied the behavior of the structure factor, changing the temperature, for / < 50. As for SP2, we 
calculated S(q) on the left and on the right side of the ideal glass region (MD simulations). 

For / = 50 the trend is the same observed for SP2, i.e. on decreasing the temperature for 77 = 0.314 the main peak 
of the structure factor decreases, showing a loss of the correlation between particles. The opposite trend is observed 
on the right side (77 = 0.628). 

For / = 46 the behavior of the structure factor is different. In Fig. 1121 -a. we present a magnification of the main peak 
in the structure factor for / = 46 and r/ = 0.55 (right side of the glass region). The picture shows a decrease of the 
main peak of S(q) lowering the temperature. Moreover in Fig. 1121 b we present a magnification of the main peak in 
the structure factor for the same arm number and r\ = 0.51, (left side of the glass region). The trend is exactly the 
same as for 77 = 0.55: we stress that, introducing attractive forces, star polymer solutions of 46 arms, show a loss of 
correlations among particles which are responsible for the peculiar MCT predictions fFig. H0|) . This effect is enhanced 
as the interaction is increased. 

Moreover we also studied the behavior of the structure factor for / = 45 and 77 = 0.51, in this case the trend is 
represented by a decrease of the main peak lowering the temperature. We observed that the height of the main peak 
introducing the attractive contribution and increasing the intensity of the attraction is in any case smaller than the 
peak height relative to V rep . We recall that for / = 45, MCT does not give a glass phase neither for V rep nor for SP3. 

Summarizing we can argue from Fig. 1101 that for the model interaction SP3 when rj < 0.5 the attraction seems 
always to determine a destabilization of the cages, while for packing fractions larger than 0.5 the attraction facilitates 
the formation of the glass for higher density with respect to V rep . For 77 ~ 0.5 (corresponding to the region of interest 
for / = 46) the effect of the attraction seems to determine a destabilization of the cages only. 

In Fig. 1131 we present the fluid-fluid phase diagram and the glass transition lines for / = 50. 

We conclude this section remarking that, so far, we discuss about the behavior of the first peak in S(q) only, while 
the second peak in the structure factor does not change appreciably in all the cases investigated. 

IV. CONCLUSIONS AND DISCUSSION 

In this work, we studied the structural arrest and the dynamics in star polymer solutions when attractive forces 
between macroparticles are present. The model potential we used to describe the interactions presents an ultrasoft 
repulsive term of cntropic origin at short range plus an attractive interaction at longer range. Due to the addition 
of the attractive contribution between star polymers the repulsive core becomes softer on lowering the temperature 
and consequently increasing the intensity of the attraction. We analyzed three different forms of the pair interaction 
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between stars, considering attractive forces of shorter and shorter range. In this paper we focused on stars with several 
arm numbers (from low value, / = 24, 32 to higher one / = 70) at several values of density and temperature. 

We examined the structure of the solution solving MHNC closure and performing extensive MD simulations, in this 
way we tested the accuracy of the MHNC to describe the properties which arise from such choice of the interaction 
model. Over the full range of densities and temperatures of interest MHNC and MD simulations are in a good 
agreement. Hence we have been able to conclude that MHNC is a good approximation to study systems described 
by ultrasoft-core repulsive interactions, with and without attractive forces between macroparticles. Having tested the 
accuracy of the MHNC approach, we focused our attention on the location of the ideal glass transition, studied within 
MCT. 

In particular for the case characterized by the longer range of the attraction, named SPi , the ideal-glass transition 
line we obtained is not significantly modified in comparison with the one obtained for the purely repulsive potential. 
Indeed both the structure and the dynamics of the solution are not significantly modified by introducing the attractive 
contribution. On the other hand the two other systems characterized by a shorter attractive range show some 
significant differences with respect to the purely repulsive pair interaction. 

F° r / ^ 50 SP2 and SP3, show the same qualitative behavior. In particular we can distinguish two regions in 
density, a low density regime (rj < 0.5) and a high density regime (77 > 0.5) where the system, in both cases, behaves 
differently. For low densities, on lowering the temperature, a liquid-stabilizing effect due to the attractive forces 
sets in, so that the liquid-glass transition line moves to larger values of the density. Indeed, this is an interesting 
effect since it presents the possibility to pass from a glass phase to a liquid phase decreasing the temperature. It is 
perhaps worth making a few remarks about this issue. It is now well established that systems characterized by a step 
repulsion and a short range attraction, possess a recntrancc in the the dynamical arrest. In these systems it is possible 
to melt the glass by lowe ring the temperature. This is a phenomenon now well established in theory, simulations and 
experiments [291 153 . l59l l60i |plT | . Clearly, one would be tempted to relate our findings to this phenomenology of 
colloidal solutions; however we showed that the origin of such effect in our case is different. For short range attractive 
colloidal systems at low temperature there is a glass phase originated solely by attraction and the reentrant melting 
arises from the competition between the high temperature and low temperature regime. For our system, however, 
there is no indication of an attraction dominated glass at low temperature and the destabilization of the glass is due 
to the destabilization of the cages that, in the high temperature regime, are responsible of the arrest. For 77 > 0.5, we 
found that the glass-fluid line moves to higher densities when the temperature decreases. As consequence increasing 
the intensity of the attraction, for a fixed density and arm number, the system moves from a liquid to a glass. This 
behavior is closely connected to the ultrasoft nature of the pair interaction. We recall that stars in a good solvent 
shows a rcmelting of the glass phases (or solid phases at the equilibrium) for very high densities. The addiction of 
attractive forces moves the glass-fluid line to higher densities. We could explain this effect as a contribution of the 
attraction to stabilize the cages and then to inhibit the remelting of the solution. 

For / = 45,46 the behaviors of SP2 and SP3 are very different. In the case / = 46 for SP2 we observed the 
qualitative behavior found for / > 50, while for SP3 there exists a value of the temperature, T* ~ 0.6, below which 
the system does not present any sign of a glass transition. For / = 45 considering SP2 we observed the presence of 
a glass transition, while for the simply repulsive interaction this transition is not present for such value of /, i.e. the 
attraction seems to facilitate the formation of the glass. In the SP3 case we did not observe glass phases for all the 
temperature values investigated. So in this case the effect of the addition of the attraction is a destabilization of the 
cages only. 

To verify the MCT results we performed the analysis of the dynamics by means of MD simulations. As introduced 
in section UlI Al the range of density which we examined is not sufficiently close to the ideal glass lines. For density 
closer to the transition lines (mainly for / = 45, 46) the system show a strong tendency to the crystallization. So we 
verify our MCT predictions for / > 50. 

From the whole analysis performed, we evidenced a value of the packing fraction, i.e. 77 = 0.5, which marks a change 
in the behavior of the solution. This value could be traced back to the cross-over between the two different functional 
forms of the repulsion: the logarithmic form and the Yukawa one. For star polymers in a good solvent 77 = 0.5, the 
so called overlap packing fraction, corresponds to the packing fraction above which the radial distribution function 
show a coordination shell inside the logarithmic core. In other words for 77 > 0.5 stars start to interpenetrate widely. 
Also in this case the cross-over designates approximatively the transition of the system through two different regimes. 
For rj < 0.5 the response of the system to the introduction of attractive forces is a destabilization of the cage. This is 
mainly due to the change in the repulsive contribution at short distances which become softer and softer considering 
respectively SPi, SP2 and SP3. For 77 > 0.5 stars interpenetrate more and more. The effect of the attraction seems to 
inhibit the remelting of the solution and the glass-fluid line move to higher densities. The shift of the transition line 
could be understood considering the effect of the attraction on the second shell of neighbors. Indeed if we look the 
radial distribution function very close to the glass-fluid transition line, at low temperature, we observe that the second 
shell of particles is around r = 2. In this region, see Fig. we can observe that the system feels stronger attractive 
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forces passing from SP\ to SP2 and then to SP3. We conjecture that this attraction on the second neighbor shell is 
at the origin of the extended stabilized region of the glass phase for higher densities. 

We have to remark some details relative to our analysis: first of all Vtot docs not present considerable repulsive 
maxima for large r when the temperature is lower than T* < 1. for SP2 and T* < 0.8 for SP3. Obviously the 
repulsive contribution at longer range is more pronounced for higher temperature (when the intensity of the attraction 
diminishes) and low arm numbers (when the Yukawa pair interaction decay to zero very slowly) . Anyway, also below 
these temperatures, where the pair interaction shows small repulsive contribution at large r (around ~ 3.), we verify 
that our MCT predictions do not depend on the existence of the repulsive shoulder. Indeed we performed the same 
MCT calculation with a truncated potential in which the repulsive shoulder is suppressed. We conclude that for the 
temperature and densities of interest in our study the presence of this small repulsive contribution at large distance 
does not alter the picture of the glass regions. We emphasize that the presence of the shoulder is not an effect of 
our particular choice of the attractive contribution: due to the ultrasoft-core repulsion and to the Yukawa repulsive 
contribution at long range the introduction of attractive forces of shorter range (i.e. depletion forces) could determine 
an additional repulsion outside the core. 

Since the potential presents an ultrasoft-core interaction, it is difficult to determine a natural scale of energy 
(and length scale). We decided to rcscale the temperature with respect to the integrated intensity of the attractive 
contribution (T new ) in such a way to compare our results in a more significant way. 

In Fig. E| we present the comparison between the MCT data obtained for 45 < / < 50 considering the SP2 and SP3 
models. As we can observe for f — 50 the effect of destabilization of the cages as well as the inhibition of the melt is 
more accentuated for SP3. 

Finally SP2 and SP3 show an apparent conflicting behavior concerning / = 45, 46. Due to the very small range of / 
it is very difficult to understand this peculiar behavior starting from the pair interaction. As we remind above, 7/ = 0.5 
is a value for which star polymers in athermal solvent, at the equilibrium, are in a solid phase and in particular 0.5 
is just in the middle of the solid region (34 < / < 70). Around this region we find, that the system has a very strong 
tendency to crystallize. Indeed, despite the change observed in the structure factor, and considering the connection 
between equilibrium phase diagram and glass transition line f|30|). we expect for SP2 and SP3 a small shift of the 
solid region with respect to the purely repulsive case. Notice that this region corresponds to the glass region we 
determined for stars with / = 45,46. It is not surprising that our data in the case of / = 45,46 show an ambiguous 
trend between SP2 and SP3. For the above mentioned reasons we do not believe it is very interesting to analyze 
further on / < 46. 

We conclude from this analysis that the details of the phase diagram, concerning glass transition, is very sensitive 
to the particular, specific form of the attractive contribution. This could be mainly connected to the change of the 
short range repulsive contribution when attractive forces between stars are introduced and not dependent on the 
choice of the specific model. We expect the same behavior also considering depletion forces. Moreover in the case of a 
general mixtures of micelles with hydrophobic group at the ends of the polymer chains, as discussed in section^J the 
attractive contribution will be at very short range and close to interparticle separation equal to a. This determines 
a change in a repulsive interaction in the region where in star polymer solution there is a cross over between the 
logarithmic form and the Yukawa one. We decide to complete our analysis considering a specific system described 
by specific attractive forces. Starting from this analysis a future perspective of our work on star polymer solutions 
is a more direct comparison between theory and experiments about the origin and the description of the attractive 
interaction. In this sense we decided to turn our attention, more in general, on systems of micelles which can be 
described by soft-core potential plus attraction at shorter range, similar to star polymer macromolccules. There 
exist many reasons to further study star polymer solutions: a precise understanding of their properties will give the 
possibility to make progress in the exploration of 'hybrid' polymeric-colloidal materials such as irregular multiarm 
stars, self-organized structures resulting from telechclic linear homopolymers and copolymers, polyclcctrolytc brushes, 
micelles with chemically fixed cores. 
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FIGURE CAPTIONS 



FIG. 1. f3Vtot(r) versus the interparticle separation r/a for two different temperature and / = 50; T* = ksT/C. 
We call SPi the interaction studied in ref.J^ i.e. A = 2.1a, B = 0.35a, SP2 the interaction corresponding to 
A = 1.95<7 and B = 0.21<r, finally SP3 corresponds to A = 1.875<r and B = 0.155<r. In figure we also show V rep 
corresponding to / = 50. 

FIG. 2. Athermal solvent: comparison between the main peak of the structure factors we obtained with MHNC 
closure and MD simulation, with the results obtained by Watzlawek et al. by means of MC simulation and RY 
closure; -q = 0.60, / = 32. (S(q) versus qa) 

FIG. 3. Comparison between the radial distribution function (top) and the structure factor (bottom) we obtained 
with the MHNC closure (solid line) and MD simulation (opaque circles); ?y = 0.628, T* = 0.25, / = 50, SP 2 . 

FIG. 4. SP\ fluid-fluid phase diagram (opaque triangles up) for a star polymer solution calculated by means of mean 
field theory and the ideal glass transition lines (filled circles) determined with MCT, which delimit the region 
where the system is a glass. The inset shows a magnification of the coexistence curve at higher density 3l| . On 
the top of the main figure we present some results about the purely repulsive interaction, which corresponds 
to the limit of very high temperature: two successive squares (from low to higher densities) delimit the regions 
where the system is solid at the equilibrium (by Watzlawek et al., MC simulation ^lj|); the filled triangles up 
represent the densities which delimit the glass region. Notice that the second critical point survives with respect 
to the glass transition while the region around the second critical point becomes metastable with respect to the 
freezing. Moreover we can observe that the densities which delimit the glass region are effectively the same with 
and without attractive interaction. 

Lines are simply a guides to the eye. 

FIG. 5. MCT fluid-glass and glass-fluid lines (computed with MHNC) for different / values and A = 1.95cr, B = 
0.21c (SP2): reduced temperature T* = ksT/C versus reduced packing fraction 77 = -|/0<7 3 . Opaque symbols, 
circles, triangles, squares, stars and diamonds, correspond respectively to / = 45,46,50,58,70. In the region 
between two line, fixed /, the system is in a glass phase. The equivalent symbols on the top of the figure delimit 
the glass regions for V re p{T). We recall that in the limit of T very large we return to the simply repulsive system. 

Lines are simply a guides to the eye. 
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FIG. 6. SP2: noncrgodicity parameter f q versus qa for stars with / = 50, 77 — 0.8. Notice the smaller width of this 
parameter for high temperatures. 

FIG. 7. SP2, / = 50: magnification of the structure factor (main peak) we obtained with MD simulations at different 
temperatures, a) r\ = 0.314, b) r\ = 0.628; / = 50. 

FIG. 8. Reduced diffusion coefficient (MD) for different / values and A = 1.95a, B = 0.21cr (SP 2 ). The symbols, 
plus, circles, triangles, squares, correspond respectively to / = 42,45,46,50 (left panel rj = 0.366, right panel 
T) = 0.628). 

Lines are simply a guides to the eye. 

FIG. 9. SP2'. fluid- fluid phase diagram for a star polymer solution calculated by means of mean field theory and the 
ideal glass transition lines determined with MCT. The legend is the same that in Fig.0J Notice that the second 
critical point survives with respect to the glass transition while the region around the second critical point 
becomes metastable with respect to the freezing. Moreover we can observe that decreasing the temperature the 
densities which delimit the glass region move to higher values, (in these cases the effect is small). 

Lines are simply a guides to the eye. 

FIG. 10. MCT fluid-glass lines (computed with MHNC) for different / values and A = 1.875a, B = 0.155er (SP 3 ): 
reduced temperature T* = ksT/C versus reduced packing fraction 77 = pa 3 . Filled triangles, squares and 
diamonds correspond respectively to f — 46, 50, 70. The equivalent symbols on the top of the figure delimit 
the glass region corresponding to V rep (r). We recall that in the limit of T very large we return to the simply 
repulsive system. 

FIG. 11. Reduced diffusion coefficient (MD) for different / values and SP3. The symbols, plus, circles, triangles, 
squares, correspond respectively to / = 42,45,46,50 (left panel 77 = 0.366, right panel r/ = 0.628). 

Lines are simply a guides to the eye. 

FIG. 12. SP3, f = 46: magnification of the structure factor (main peak) we obtained with MD simulations at 
different temperatures, a) 77 = 0.55, b) 77 = 0.51. 

FIG. 13. 5*^3: fluid-fluid phase diagram for a star polymer solution calculated by means of mean field theory and 
the ideal glass transition lines determined with MCT. The legend is the same which in Fig. 01 
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Notice that the second critical point survives with respect to the glass transition while the region around the 
second critical point becomes mctastable with respect to the freezing. Moreover we can observe that decreasing 
the temperature the densities which delimit the glass region move progressively to higher values. 

Lines are simply a guides to the eye. 

FIG. 14. Magnification of the low temperature MCT fluid-glass lines (computed with MHNC) for different / values, 
SP2, SP3. The temperatures are rescaled with respect to the integrated intensity of the attractive contribution. 

Lines are simply a guides to the eye. 
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